A LOW-REYNOLDS-NUMBER TREADMILLING 
SWIMMER NEAR A SEMI-INFINITE WALL 

KIORI OBUSE* AND JEAN-LUC THIFFEAULTt 



Abstract. We investigate the behavior of a treadmilling microswimmer in a two- 
dimensional unbounded domain with a semi-infinite no-shp wall. The wall can also be 
regarded as a probe or pipette inserted into the fiow. We solve the governing evolution 
equations in an analytical form and numerically calculate trajectories of the swimmer 
for several different initial positions and orientations. We then compute the probability 
that the treadmilling organism can escape the vicinity of the wall. We find that many 
trajectories in a 'wedge' around the wall are likely to escape. This suggests that inserting 
a probe or pipette in a suspension of organism may push away treadmilling swimmers. 



1. Introduction. The locomotion of microorganisms is an active re- 
search area of fluid dynamics and biology (see for instance the reviews 



10 15]). As their motion occurs on very small length scales and speeds, 
their dynamics is governed by low-Reynolds-number hydrodynamics, where 
inertial forces are negligible in comparison to the viscous effects of the fluid 
(Stokes flow). 

Many studies deal with such dynamics in unbounded or very large do- 
mains [T T2jT8]. In reality, however, most organisms are in the vicinity 
of other bodies or boundaries, where hydrodynamic interactions can have 
a significant effect on their motion. The importance of boundaries has 
also been suggested by experimental observations. For example, some re- 
search suggests that microorganisms are attracted to solid walls [3 17 19 



Berke et al. [2] measured the steady-state distribution of E. Coli bacte- 
ria swimming between two glass plates and found a strong increase of the 
cell concentration at the boundaries. They also theoretically demonstrated 
that hydrodynamic interactions of swimming cells with solid surfaces lead 
to their reorientation in the direction parallel to the surfaces. Lauga et 
al. [9] showed that circular trajectories are natural consequences of force- 
free and torque-free swimming and hydrodynamic interactions with the 
boundary. This leads to a hydrodynamic trapping of the cells close to the 
surface. Drescher et al. [6] found that when two nearby Volvox colonies 
swim close to a solid surface, they attract one another and can form stable 
bound states in which they 'waltz' or 'minuet' around each other. These ob- 
servations suggest that, in order to obtain a comprehensive understanding 
of low- Reynolds-number locomotion, it is necessary to study hydrodynamic 
interactions between microorganisms and boundaries. 

Some of the phenomena stated above have already been confirmed 
by numerical simulations |8l|16|. However, not many physical explana- 
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tions have been given to the locomotion of microorganism near boundaries. 
Berke et al. [2] have captured the swimming microorganisms' attraction to 
boundaries by modehng the swimmer as a force dipole singularity. However, 
contrary to the experimental findings, the microorganism in this model 



crashes into the boundary in finite time. Or and Murray 14 studied the 
dynamics of low-Reynolds-number swimming organism near a plane wall. 
They analyzed the motion of a swimmer consisting of two rotating spheres 
connected by a thin rod as a simple theoretical model of a 'treadmilling' 
swimming organism. They found that when the spheres rotate with dif- 
ferent velocities their model has a solution with steady translation parallel 
to the wall, and that under small perturbation the swimmer exhibits a 
'bouncing' motions parallel to the wall. These results have recently been 
verified experimentally on a macroscale robotic prototype swimming in a 
highly viscous fluid 20 . Furthermore, Crowdy and Or |4] have proposed 
a singularity model for swimming microorganisms placed near an infinite 
no-slip boundary. Their model was based on a circular treadmilling organ- 
ism which has no means of self-propulsion, that is, the organism doesn't 
move unless it interacts with a boundary. (For example, the organism 
may be creating a feeding current.) They proposed an appropriate Stokes 
singularities that represent the flow field created by this treadmilling or- 
ganism. By studying the interaction between these singularities and the 
no-slip wall, they formulated explicit evolution equations for the motion 
of the organism, and fully characterized its motion near the wall. They 
found trajectories with a periodic bouncing motion along the wall which 



had remarkable similarity to the trajectories shown in Or and Murray 14 
Crowdy and Samson [s], using the point-singularity model, investigated the 
dynamics of treadmilling organism near an infinite no-slip boundary with a 
gap of a fixed size. They employed a conformal mapping technique to avoid 
the difficulty in treating the image of the treadmilling organism on the wall. 
They found that the treadmilling organism can exhibit several qualitatively 
different types of trajectories: jumping over the gap, rebounding from the 
gap, being trapped near the gap, and escaping the gap region even when the 
organism has initial position in the gap. They also performed a bifurcation 
analysis in terms of the model parameters, and demonstrated the presence 
of stable equilibrium points in the gap region as well as Hopf bifurcations 
to periodic bound states. This reduced model also exhibited a global glu- 
ing bifurcation in which two symmetric periodic orbits merge at a saddle 
point into symmetric bound states having more complex spatio-temporal 
structure. 

In the present paper we examine the dynamics of a treadmilling or- 
ganism near a semi-infinite no-slip wall, modeled as a flat plate of zero 
thickness. Though this is a special case of the model of Crowdy and Sam- 
son [5] , it deserves separate investigation because of the simpler equations 
involved, and because the semi-infinite wall can be regarded as a probe or 
pipette inserted in the system, a common situation in microbiology. We 
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also analyze the trajectories in a very different manner to [s], as we attempt 
to quantify the probability of escape from the wall's vicinity, assuming the 
treadmillers are randomly oriented. 

2. Model of a treadmilling microorganism. Following previous 
authors [4j|5] , we consider a two-dimensional model for a microorganism in 
the {x, y)-plane, which we treat as the complex plane with z = x + iy. Our 
derivation is a special case of [S] , who considered a microswimmer near a slit 
or gap in an infinite wall. Nevertheless, as mentioned in the introduction, 
the semi-infinite wall is important enough to be treated separately, as it 
arises in the neighborhood of a probe or a pipette. 

The Stokes equations which describe the motion of an incompressible 
viscous fluid are 

Vp = TjAu, V • M = 0, (2.1) 

where A is the Laplace operator, u = {ux,Uy) is the fluid velocity, and 
p and T] are the pressure and dynamic viscosity, respectively. As we are 
considering a two-dimensional flow, we can introduce a stream function 
such that the velocity is given by = dip/dy, Uy = —dip/dx. Then the 



Stokes equations (2.11 reduce to the biharmonic equation A^V^ = 0. The 
complex velocity is W = Ux + iuy = —2idip/dz, with 



W^ux+iuy = -2i-^ = f{z) + zf'iz) + (2.2) 
^ = lu,[zf{z)+g{z)]. (2.3) 

where f{z) and g{z) are called Goursat functions; they are analytic every- 
where in the flow domain, except where isolated singularities are introduced 
to model swimmers. For a treadmilling swimmer, we take f{z) to have a 
simple pole at z = z^;, so that 

f{z) = ^^ + ..., g'(z) = ^J^ + .... (2.4) 
where /x € C and the form of g' is forced by the requirement that the 



complex velocity (2.2) have no higher than a |z — z^j ^ singularity. This 



solution corresponds to a stresslet of strength n at z^. The ellipses in (2.4 1 



indicate analytic terms. The expansion (2.4) is the basic solution for a 
treadmilling swimmer, which does not have any self-propulsion in itself, 
but can move due to its interaction with boundaries |4, 5, 11, 14j[20j. 

In a simple model, we assume that the treadmilling organism has a 
circular body of radius e, with a center at Zd{t) = Xd{t) + iyd{t)- We also 
assume that, with respect to the angle 9{t) of the head of the treadmilling 
organism from the real axis, surface actuators of the treadmilling organism 
induce a tangential velocity profile given by ,4,| 



U{(j),e) = 2Vsm{2{(j)- 9)), 



(2.5) 
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where ^ is a constant and (j) is the angle measured from the positive x 
direction. 

Next consider the treadmiller near an infinite wall along the x axis. To 
satisfy the no-slip boundary condition at the wall, we make the expansion 

fiz,t)^ ^ + f^ + {z-zd{t))h + -- - , 

a (2.6) 

f{z,t) having no Stokeslet term and g{z,t) having no rotlet term implying 
that the treadmilling organism is force-free and torque-free. We use the 
boundary condition to find 

/.i = —zee, a = ^Zd, b = fie^ — ice^ — (2-7) 

where c{t) = —iVcxp{—2i6{t)). We set the time scale of the motion by 
letting V = e^^ so that /i(t) — exp{2i9{t)). The coefficients /o, /i, and go 
are given in |4]. 

The time derivative of the position and orientation of the treadmilling 
organism is obtained by equating the time rate of change of position to the 
local fluid velocity, and the time rate of change of orientation to half the 
local vorticity: 

^^-.fo + zji+go, f = -2Im/i. (2.8) 



Equation (2.8 1 can then be solved as a set of three ODEs determining the 
motion of the treadmiller. 

Now we turn to a semi-infinite wall, extending along the negative x 
axis. The conformal mapping ^ = iz^^'^ maps the z plane to upper-half C 
plane, with the negative cc-axis of the z plane mapped to the real axis in 
the plane. In the ^ plane we can use a similar singular expansion as ( |2.6[ ), 
but we must take care to map the boundary conditions to the C plane. We 
omit the lengthy details, which are similar to Crowdy and Samson [sj. 



See 13 for a more complete derivation. All that is required for simulating 



the swimmer trajectories are the coefficients that appear in (2.8), which 
are 

e^fi (2|zd|2 - 2z2 - 3e2) ^ {2\za\^ + 2zj - fl 
Jo ~ ' ' 



4Zrf 423^3/2 8Z2z2 g^Z^^ 



(2.9a) 



1 / 3a. 9z'J\'fl 3zT{2\z,\'-2zj-3e^)fl 

-izT (2|zd|2 + 2z2-3e2)/2^ 



5/2 



(2.9b) 
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(a) Trajectory: 

(xo.yo. 0o) = (-10, 10, 37r/4) 



(b) Trajectory: 

(xo.yo, 00) = (-10, 10, 187r/25) 

I5| ' ' ' ' 




Fig. 1. Examples of trajectories from the initial points (aj^joij/do) marked with 
circles, with an initial orientation Bq. A treadmilling organism can escape from the 
wall, such as in (a) and (d), end up above the wall as in (b), or underneath the wall as 
in (c). 



and 

3^Zrf lOe^// 3e^/i [zd - Zd) Jl 



90 



16zj 32z3 SZ^Zd 4^3^1/2 

(2|zrf|2 - 6zj ~ 3e2) ft 



16Z2z 



.5/2 



(z-y^ + z), (2.9c) 



where Z = z^^ + z^"^ . 



3. Results of numerical simulations. We now present the results 



of numerical simulations of the governing evolution equations (2.8), to- 



gether with the coefficients (2.9). The radius of the circular treadmillers 
is set to e = 1, giving the reference length scale. The time scale is set 
hyV^ £-1 in 

Figure [T] shows examples of swimmer trajectories for different initial 
conditions for position Zdo and angle ^o- Some trajectories, such as (a) and 
(d), end up far away from the wall; we refer to those as escaping trajectories. 
Others, such as (b) and (c), remain close to the wall for all time, exhibiting 
the 'bouncing' behavior noted in 4 . A fourth type of trajectory (not 
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Fig. 2. (a) The x-y plane divided into five regions. The thick black line corresponds 
to the wall, (b) An example of a pie chart. The different colors as a function of initial 
orientation 0o indicate which region of (a) the treadmilling organism ends up in after 
a sufficiently large time. 

shown) crashes into the wall, but this is an unphysical consequence of the 
boundary condition at the swimmer's surface only being approximately 
satisfied. It has been verified numerically that the qualitative features of 
the trajectories are not affected by this approximation |7]. 

Both experimental observations and previous theoretical studies sug- 
gest that, when there is a no-slip wall near a treadmilling organism, the 
organism tends to be attracted to its own image and move towards the 
wall [2||4[[5{[8{[9| [T6{[T7|[T9] . This behavior is clearly seen at an early stage in 
all the trajectories in Fig. [l] Nevertheless, in the cases shown in Figs, [ija) 
and (d), the treadmilling organism moves away from the wall after it has 
come close to the edge of the wall. 

We now focus on the large-time behavior of treadmillers. To study 
this, we first introduce five regions shown in Fig. [2j[a). Regions 1-4 are the 
usual quadrants of the plane, but with a neighborhood of size 0.2e around 
the wall removed; this removed neighborhood is region 5. The regions 
1-5 are colored red, blue, orange, green, and white, respectively. From 
these regions, we can make a 'pie chart' around each point in the plane, 
an example of which is shown in Fig. [2j^b) . The sectors of the pie chart 
correspond to a range initial angle Oq; the color of the sector describes which 
region the treadmiller ends up in for large times (here t = 1500). Region 
5 corresponds to 'crashing' trajectories. Thus, for a given position, the 
treadmiller may end up in different regions, and the relationship between 
angle and final region is not simple. 

Figure |3] shows pie charts for many initial points {xdo,ydo)- The most 
notable feature is the complexity of the pie charts for initial points near 
the wall and with large negative x coordinates. This reflects the fact that 
the treadmilling organism changes its heading direction and the quality of 
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Fig. 3. Pie charts at t = 1500 for several initial conditions. Each pie chart is 
centered on the initial condition it corresponds to. The thick black line represents the 
wall. Note that region 3 only appears as very thin slivers in the pie charts along the x 
axis. 



its trajectory significantly when it has come to the vicinity of the edge of 
the wall, X = y = Q (depending upon its (x^, yd, 0) at the time, and so even 
a very small difference in initial condition can result in a huge difference 
to its trajectory. (This is an example of chaotic scattering.) However, pie 
charts tend to become rather simple for larger y^o for f^ny fixed x^o, since 
the treadmiller is then less influenced by the wall. 

Now let us consider the escape probability^ Psixdo, Udo), the probability 
that the treadmilling organism can escape from the wall region. We define 
this as the probability that the treadmilling organism be in region 1 or 4 
at sufficiently large time. We assume that the initial angle is uniformly 
distributed. The escape probability PE{xdo,ydo) then corresponds to the 
fraction of angles in a pie chart that are red or green, for a given initial 
position (xdOjUdo)- The escape probability for t — 1500 for different initial 
conditions is shown in Fig.|4] where each {x, y) coordinate corresponds to an 
initial position of the treadmiller. Note that for almost all initial conditions 
there is some finite probability of escaping or being trapped, though the 
escaping probability approaches unity along the positive x axis, and goes 
to zero for large negative x and large y. Observe the strong 'wedge' of 
trajectories that are likely to escape {Pe > 0.6), though there is also a 
backwards-facing wedge of trajectories that have a reasonable change of 
escaping {Pe > 0.3). This suggests that a probe or pipette inserted in a 
medium is likely to 'push away' treadmillers to some degree. Note also the 
'tongues' of abnormally-high {Pe — 0.4) escape probability near the wall. 
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Fig. 4. Escape probability at t = 1500 for different initial points. Each point 
corresponds to a given initial condition (x^jo, j/^q); integrated over all possible starting 
angles Oq . A trajectory escapes if it ends up in regions 1 or 4, as defined in Fig. [2| The 
excluded white area near the wall corresponds to initial conditions where the swimmer 
would be partially inside the wall. 



which reflects the complicated structure of the pie charts in that region. 

4. Discussion and conclusions. Following [4][5] , we derived the gov- 
erning evolution equations for a treadmilling microorganism near a semi- 
infinite no-slip wall. This can also be regarded as a two-dimensional model 
of a probe or pipette inserted in the fluid. We then numerically calculated 
the trajectories of the organism for different initial conditions. The tread- 
milling organism was usually attracted to the wall for early times, but for 
later times often escaped the wall region. Typically this happened when 
the organism came close to the tip of the wall. To investigate this behavior 
further, wc looked at 'pie chart' diagrams, where the sectors indicate which 
region a swimmer eventually ends up in as a function of its initial angular 
orientation. Initial points with larger negative x coordinates have more 
complex pie charts compared to those with smaller negative x coordinates, 
because they are sensitive to 'chaotic scattering' off the tip of the wall. 

We then examined the escape probability Pe , the probability that the 
treadmilling organism can escape to the right of the wall region, assuming 
that it is initially randomly oriented. There is an evident 'wedge' of initial 
conditions that is likely to escape the wall. This suggests that a probe or 
pipette inserted could 'push' the organisms out of the way (with the wedge 
replaced by a cone in three dimensions), though since the organisms slow 
down as they get further away from the wall the effect might not be very 
pronounced. 
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